Addition energies of a Quantum Dot with harmonic electron-electron interactions 
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We study a two dimensional electron system in a parabolic confining potential and constant magnetic 
field for the case of harmonic electron-electron interaction. We present analytic results for the 
electrochemical potential versus magnetic field and discuss the effects of correlation in connection 
with the addition energy of a Quantum Dot with few electrons. 
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Recent achievements of nanolithography in semicon- 
ductor technology allow for the fabrication of devices in 
which a definite number of electrons are confined within 
two-dimensional islands of size as small as tenths of 
nanometers [Quantum Dots (QD)] [Q. In the last few 
years there has been growing interest in the study of these 
devices in view of improving our understanding of corre- 
lated electron systems. In fact, QD are unique with re- 
spect to other structures, e.g., macromolecules and clus- 
ters, because a dot can be connected to sources and/or 
a measuring apparatus via contacts. This possibility al- 
lows investigation of the system with probes changing the 
number of particles j^]. Indeed, quite recently Tarucha 
et al. |s| have measured the tunneling current in gated 
vertical quantum dots as a function of a magnetic field 
B applied parallel to the current. In the Coulomb block- 
ade regime and in presence of a very small voltage bias, 
the current shows a sequence of peaks that occur when- 
ever the gate voltage V g is proportional (via a voltage-to- 
energy conversion coefficient) to the chemical potential 
Hn = E(N) — E(N — 1) for adding one more particle to 
the dot. Here E(N) is the ground-state (GS) energy of 
the dot once N electrons are localized in it. 

General features of the current peaks are qualitatively 
reproduced by assuming that electrons are confined by 
a parabolic potential of frequency luo and by adding 
the charging energy E c h — VqN(N — l)/2 to account 
for the electron-electron repulsion [Constant Interaction 
(CI) model] H,|[. The GS energy is obtained, for any 
value of B, by filling the lowest one-particle free har- 
monic oscillator levels with electrons of both spin. In 
this model the observed increasing of the addition en- 
ergy A N = E(N + 1) + E(N — I) — 2E(N) for N = N p , 
where N p — 2,6,10,20,..., is easily related to the shell 
structure of the two dimensional harmonic potential spec- 
trum, which leads to a marked increase of [In whenever 
a new shell is opened. The CI model then reproduces the 
oscillations of the current peaks with the field B as well 
as their shift in pairs. Because the CI GS wave function 
has minimum total spin, its agreement with the experi- 
mental pattern also confirms that in some QD, like the 
one in Ref . || , the effective g-factor can be very small 
||. However, there are some evident features of the ex- 
periments that cannot be reproduced by the simple inde- 



pendent electron CI Hamiltonian. Because the harmonic 
single-particle levels are equally spaced, according to the 
CI model (at B = 0) An should be constant (An = Vq) 
within a unfilled shell. For the same reason one should 
then find An — Vq + huiQ for any N p . While the exper- 
imental pattern in Ref. fl allows a clear identification 
of a shell structure for dots of less than w 20 electrons, 
it also clearly shows a smooth shrinkage of the spacing 
Ajv within a unfilled shell with growing N, as it was 
previously noticed in single electron capacitance spec- 
troscopy H . Then, one observes that the expected peaks 
at An p decrease in height each time a new shell opens. 
In general, as discussed by Schmidt et al. 0, fittings of 
capacitance data within the CI model indicate that the 
confinement strength of the potential appears to decrease 
rapidly with increasing energy. Question arises whether 
such an effect can be attributed to electron-electron cor- 
relation. 

In this letter we consider the model of harmonic in- 
teractions (HI) between TV electrons S in a parabolic 
confining potential. Although there is tenuous justifica- 
tion for this model interaction, when the magnetic field 
is not too small and the electron number not too large we 
show that it embodies correlation effects which correctly 
reproduce interesting experimental features. The model 
is well known and its popularity relies on the fact that 
for fully spin polarized electrons the Laughlin trial wave 
function, which is successful in describing the Coulomb 
gas in a quantizing magnetic field, has a form similar to 
one of its eigenstates ||. Successively, the model has 
been discussed, again in the fully polarized spin sector, 
in connection with the QD problem in a large magnetic 
field orthogonal to the dot plane |l(| . In view of the find- 
ing w 0, results about the maximum spin sector are 
not relevant to the present discussion. We present some 
details of the solution because little information can be 
found in literature about the GS wave function in an ar- 
bitrary spin sector. 

We do not expect the HI model useful in discussing 
the region of small magnetic field. In fact, there are evi- 
dences H that in this case the spin sector of the QD GS 
obeys Hund's rule, whereas the GS of the HI model at 
B = has necessarily lowest spin Instead, for higher 
magnetic field our results closely reproduce features of 
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the experimental data of Ref. |3|] and are summarized in 
Fig. (1). The narrowing, compared with the CI model, of 
the distance between current peaks is transparent. The 
shift, at fixed Vo; of the orbital angular momentum tran- 
sitions which take place with increasing B to lower values 
of the magnetic field also clearly appears in our results. 
Then, close to the point where the oscillations in the peak 
location drop, for large enough AT, the peaks making up a 
pair have an intriguing "out of phase" behaviour, similar 
to the experimental pattern. 

The HI Hamiltonian in a constant magnetic field 
B is H B = H + hui c (L 3 + g*S 3 )/2, where HL 3 = 

J2iLi( x i,iP2,i - %2,iPi,i) and S 3 = £) i=1 S3,i are the 
third component of the total angular momentum and 
total spin operators, respectively, r.; — {x\^,X2,i) and 
p t = —ihd/dri are the coordinates and momentum of 
the z-th electron, uj c = eB/(m*c) is the cyclotron fre- 
quency, and 

*2 w „ » 2 N 

+ E W-f \ri-rj\ 2 ). (1) 

i<»<i<w 

Here m*, g„ denote effective parameters and for the 
time being we set = e = h = 1. The frequency 
lo 2 = uo 2 + uo 2 /A enters Ho upon choosing the gauge 
A = (B/2)(—X2,xi) for the vector potential and rep- 
resents the frequency of the effective parabolic confining 
potential, whereas U > is the strength of the inter- 
action. For U — then Hb reduces to the CI model. 
Introducing Ay = Q 2 5ij + UJij, where fi 2 = uj 2 — NU 
and Jij denotes the matrix with all unit entries, (i.e., 
Jij = 1 Vi,j), the potential energy entering Eq. (|l|) 

can be compactly written as V — Y^ij=i^-ij r i ' r j/2> 
so that Hq is bound from below if Ay is positive definite. 
The symmetric matrix Ay is diagonalized by any uni- 
tary matrix Ui V satisfying to Yli=i = 0, W ^ N, and 
UiN = , Vi. Its eigenvalues Ai = . . . = Ajv_i = ft 2 , 
\n = lu 2 are readily evaluated, so that positivity is en- 
sured whenever lo 2 > NU. Hence, if we limit the discus- 
sion to N not too large, the unphysical feature of deal- 
ing with an interaction unbounded at large distances is 
compensated by the presence of the confining potential. 
We then introduce normal coordinates y v — Yli=i U^ri. 
The Laplacian A = Yli=i^/^ r i ' 9/dri, the angular 
momentum L 3 , and the operator R 2 = Yli=i \ r i\ 2 are 
invariant under unitary transformations and from the 
invariance of R 2 one gets the simple but key identity 

J2u=i \Vu\ 2 = J2i<j \ r i - r j\ 2 / N - In the new basis the 
equation H B ^ = is immediately solved, because 
the key identity allows to write H as a sum of sepa- 
rated harmonic oscillators of frequencies \f\^. However, 
the straightforward normal-mode approach is of no much 



practical help, because the main problem one faces is to 
account for the identity of the particles. Eigenfunctions 
of Hb with definite symmetry under particle permuta- 
tions can be factorized in the form \P = ^ cnl ^ r , in the 
usual way JO], where $ cm is the completely symmet- 
ric harmonic wave function of the center of mass (cm.) 
coordinate r = N^^y N , with energy E cra and angular 
momentum L§ m , while the relative motion wave function 
^ r must be solution of the equation 

- + ^ J2 \ri - rj\H r + y (L 3 + g*S 3 )V r 

i<j 

= {E — E ch - E CU1 - Y L T)^r, (2) 

of same symmetry as ^ and of zero total linear momen- 
tum p = EiLiPi- 
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FIG. 1 Chemical potential fiN versus magnetic field B for 
TV < 22: (a) HI model for hU/(m*u ) = 0.07 meV; (6) CI 
model. Parameter values are %uio = 3 meV, Tiw c /B — 1.63 
meV/T, V = 1.5 meV, g, = -0.03. 



Henceforth we need only consider the lowest energy 
cm. wave function and for the time being we thus 
set "Pcj = exp{— wY^ij=i r i ' r j/(2-W)}j f° r which one 
has E'en = k> and L 3 m = 0. If one looks for eigen- 
states of Eq. (§) in the form * r = $* 0j where * = 
exp{— QJ2i<j \ r i ~ r j| 2 /(2A r )}, and employes holomor- 
phic coordinates Zj, = xu + ix%,i, and Zi = x\^ — 1x2,1, 
one gets that the unknown function <!> must satisfy the 
couple of equations 

N 
i=l 

= (E — £0 — E ch )$, P$ = 0, (3o,ft) 

where £q — Q(N — 1) + lo is the zero-point energy, 
Cl± = Q,±uj c /2, and we have used the shorthand notation 
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di = d/dzi, and d% = d/dzi. Because the function "Jo is 
completely symmetric, \l/ r and $ must have same sym- 
metry. Eq. (3a) is a sum of N separated Hamiltonians 
and its solutions are built up in terms of single-particle 
orbital wave functions /„ lir!2 , with eigenvalues £ nij n 2 , 
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(4) 



where ni and n 2 are non negative integers. Appar- 
ently, the original interacting problem reduces to a free 
problem and the main effect of the interaction seems 
just related to a redefinition of the effective frequency 
cu — > O. However, Eq. (3b) is highly non trivial and 
already the two-particle problem shows that it can be 
satisfied in general only by taking linear combinations 
of degenerate A-particle solutions. Moreover, the fre- 
quency £1 = \/lu 2 — NU depends on N, so that energy 
differences like E(N) — E(N — 1) cannot be analyzed in 
terms of one-particle levels. In these respects the system 
is correlated. 

Fortunately, as expected on general ground, one can 
easily check that the ordinary Slater determinant Z ob- 
tained by filling the N lowest states @ (so that S3 — 
N/2) is a solution of both Eqs. ( 3a, 6 ). Fig. (2) shows 
some typical shape of the set of occupied levels. By in- 
creasing B one observes a depletion of the levels in the 
SE side and an extra filling in the NW side, because 
the contribution Hl — u c L z /2 favours decrease of the 
orbital angular momentum. Noticing that the functions 
fni,n 2 are polynomials and using standard properties of 
determinants under column addition and multiplication, 
it is easy to see that for any set of the form depicted in 
Fig. (2) the Slater determinant can be rearranged into the 
compact form Z = det^" 1 ^™ 2 ]. This expression is invari- 
ant under translation z% — > Zi + zo, Vi, so that PZ = as 
requested. Introducing the spin index a = 1, 2 for up and 
down electrons, respectively, it follows that the GS wave 
function ^gs of Hb in the spin sector S3 = (Mi — M 2 )/2, 
N = Mi + M 2 , is obtained, up to a normalization con- 
stant, by antisymmetrizing the product of the spin-up 
and spin-down wave functions. Denoting A the antisym- 
metrization, we formally have 



h:\xc 



(5) 



where Zi a , with i± = 1, . . . , Mi and i 2 = Mi + 1, . . . , N, 
are the coordinates of spin- up and spin-down electrons, 
respectively, \ii Xi-, are the totally symmetric spin func- 
tions of spin-up and spin-down electrons, and (nj, ) are 
the labels of the (lowest energy) levels occupied by spin-cr 
electrons. The matrix fy = QSij + N^ 1 ^ — il)Jij is the 



square root of Ay and enters Eq. (||) upon combination 
of ^0 with the cm. GS wave function. 

For zero magnetic field the GS is obtained by filling 
consecutively the levels along successive diagonals n x + 
n 2 = k, k = 0, 1, ... with two electrons of opposite spin, 
so that ^gs has S z — or S z = ±1/2. Completion of a 
diagonal corresponds to one more shell filled. When the 
magnetic field is turned on, for a wide range of B the GS 
remains in the lowest spin sector due to the smallness of 
<7* , whereas as previously seen the shape of the Fermi sea 
modifies in order to lower the orbital angular momentum. 




FIG. 2 Occupied single-particle quantum numbers (121,712) 
of 15 spinless fermions for six different values of the magnetic 
field B (T). Parameter values are the same as in Fig. (la). 



The lowest attainable value of the orbital angular 
momentum Ef ln — — M a (M a — l)/2 is obtained 
by filling the levels (n^n^) = (0,TOo-), where m a — 
0,1,..., M a — 1. In this case the Vandermonde deter- 
minant Z = detfzf 1 ] = rii<j (% — %i) allows a closed ex- 
pression for the GS wave function ^qs- Denoting {r^, e^} 
orbital and spin coordinates of the z-th electron, Eq. (|5|) 
gives the Laughlin-like |T^] state 



1 N 

=exp{--^ ZiTijZj} 



Kj 



(6) 



Our results for the chemical potential vs. magnetic field 
are shown in Fig. (1). The curves are splitted by the 
charging energy Vq but the A-dependence of f2 leads 
to a sizeable reduction of the energy scale of the HI 
model compared with the CI model. This behaviour is 
in agreement with the experimental pattern of Ref. j^] 
and with the discussion in Ref. 0. Although one can- 
not overestimate the model nature of the harmonic in- 
teraction, our results indicate that the rapid decrease of 
the confinement strength with increasing electron num- 
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ber can be due to pure correlation effects. The curves 
(for N > 5) perform oscillations that signal transitions 
to states of lower L3 , and drop at the field value where 
the minimum angular momentum state becomes en- 
ergetically favourable. It is easy to see that this happens 
when 2Q/uj c < M 1 /{M 1 - 2) (with M x > M 2 ). In the 
CI model the curves are perfectly paired because of the 
spin degeneracy of the one-particle levels. For example, 
the transition to the state of lowest L3 takes place when 
2ui/lu c — Mi I (Mi — 2) . This relation gives the same value 
of the magnetic field at which the transition takes place 
both for N = 2M and for N = 2M-1. Instead, in the HI 
model the iV-dependence of Q breaks this perfect pairing, 
in particular for large N where becomes small. While 
usually the angular momentum of the GS is a decreasing 
function of the number of particles, at the field values 
where the paired curves come out of phase one instead 
finds L Z (2N - 1) < L Z {2N). For instance, at B = 3.6 T 
one has L z (21) = —89, whereas L z {22) = —88. Hence, 
the strong reduction of the escillator frequency at large N 
allows for level crossing even within a single pair of par- 
ticles which are added with opposite spin. This in turn 
leads to the oscillations out of phase of the peaks form- 
ing a pair, a feature which is absent in the CI model. We 
conclude that, in spite of the simplicity of the model for 
including electron-electron correlations, many features of 
the experiment of Ref. |3| are nicely reproduced. 
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